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Complete Tenth-Order QED Contribution to the Muon g — 2 
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We report the result of our calculation of the complete tenth-order QED terms of the muon 
g — 2. Our result is a^^''^ — 753.29 (1.04) in units of (q/tt)^, which is about 4.5 s.d. larger than 
the leading- logarithmic estimate 663 (20). We also improved the precision of the eighth-order QED 
term of a^, obtaining a^**' = 130.8794 (63) in units of (q/tt)*. The new QED contribution is 
a^(QED) — 116 584 718 951 (80) x 10~^*, which does not resolve the existing discrepancy between 
the standard-model prediction and measurement of a^. 

PACS numbers: 13.40.Em,14.60.Ef,12.20.Ds 



The anomalous magnetic moment of the muon has 
been studied extensively both experimentally and the- 
oretically since it provides one of the promising paths 
in exploring possible new physics beyond the standard 
model. For this purpose it is crucial to know the predic- 
tion of the standard model as precisely as possible. 

On the experimental side the current world average of 
the measured is [ll, 0] : 

a^(exp) = 116 592 089 (63) X 10"" [O.Sppm]. (1) 

New experiments designed to improve the precision fur- 
ther arc being prepared at Fermilab Q and J-PARC [1]. 

In the standard model, can be divided into electro- 
magnetic, hadronic, and electroweak contributions 

0-^1 = a^(QED) + (hadronic) + (electroweak). (2) 

At present (hadronic) is the largest source of theoreti- 
cal uncertainty. The uncertainty comes mostly from the 
0{a^) hadronic vacuum-polarization (v. p.) term, a be- 
ing the fine-structure constant. The lattice QCD simula- 
tions have attempted to evaluate this contribution 
At present, most accurate evaluations must rely on the 
experimental information. Three types of measurements 
are available for this purpose: (1) e+e" — >■ hadrons, (2) 
I' + TT^ , (3) e+e" — 5- 7 -I- hadrons. These pro- 
cesses have been investigated intensely by many groups 
UHlil.We list here one of them [l3|: 



,(had. v.p.) 6949.1 (37.2)exp(21.0),.ad x 10 
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(3) 



which overlaps other values based on the e+e" data [ll|, 
12l | and makes the standard-model prediction closest to 
the experiment (H]). The next-to-lcading-order (NLO) 
hadronic vacuum-polarization contribution is also known 



,(NLO had. v.p. 



-98.4 (0.6)exp(0.4)rad X 10- 
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The hadronic light-by-light scattering contribution (l-l) 
is of similar size as a^(NLO had, v.p.), but has a much 
larger theoretical uncertainty [lj-[l7| 



a^(had. l-l) = 116 (40) x 10" 



(5) 



where the uncertainty 40 x 10"^^ covers almost all values 
obtained in different publications. 

The electroweak contribution has been calculated up 
to 2-loop order [li-l2l|: 



a^(weak) = 154 (2) x 10 



'11 



(6) 



Since this uncertainty is 30 times smaller than the ex- 
perimental precision of ([T|), it can be regarded as known 
precisely. 

The primary purpose of this letter is to report the com- 
plete numerical evaluation of all tenth-order QED contri- 
bution to . It leads to a sizable reduction of the uncer- 
tainty of the previous estimate by the leading-log approx- 
imations [22, [2^ . We have also improved the numerical 
precision of the eighth-order QED contribution including 
the newly evaluated tau-lepton contribution. Together 
they represent a significant reduction in the theoretical 
uncertainty of the QED part of a^. 

The QED contribution to can be evaluated by the 
perturbative expansion in a/n: 



a,(QED) = ^ 



(7) 



where a^^"'' is finite thanks to the renormalizability of 
QED and can be written as 



,(2n) 



-t- A2^"-'(m^/?7ie) + A':^"' {mfj_/mr) 



(2n) 



-I- A^^"-'(TO^/771e,m^/m^). 



(8) 



,(2n) 



is independent of mass and universal for all lep- 



(4) tons. A[^\ A[^^ and Af are known exactly (2j-l2 



(6) 



p;^ 



FIG. 1. Vertex diagrams representing 13 gauge-invariant 
subsets contributing to the lepton g — 2 at the eighth-order. 
Solid and wavy lines represent lepton and photon lines, re- 
spectively. 



,(2n) 



and 



Mass dependence is known analytically for 

) fr.^ „ - raJooll We reevaluated them us- 



A^^"-' for n = 2,3 |28|-l32| 



ing the latest values of the muon-electron mass ratio 
m^/mf. = 206.768 2843 (52) and/or the muon-tau mass 
ratio mf,/mr = 5.946 49 (54) x IQ-^ In the same 

order of terms as shown on the right- hand-side of ([8]), the 
results are summarized as follows: 



.(4) 



0.5, 

-0.328 478 965 579 



1.094 258 312 (83) 



+ 0.780 79 (15) X 10"^ 
0.765 857 425 (17) , 
1.181 241 456 . . . + 22.868 380 04 (23) 
+ 0.360 70 (13) X 10"'"' -I- 0.527 76 (11) x 10"^ 
= 24.050 509 96 (32) . (9) 

The value of a/i has been obtained mostly by nu- 



merical integration |34|-l36l| . They arise from 13 gauge- 
invariant sets whose representative diagrams are shown 
in Fig. [1] We have reevaluated some of them for further 
check and improvement of numerical precision. The re- 
sults for the mass-dependent terms are summarized in 
Table H 

From the data listed in Table |T] and the value of A^*' 
from Refs. 35l437l|. we obtain the following value for the 



TABLE I. The eighth-order mass-dependent QED contribu- 
tion from 12 gauge-invariant groups to muon g — 2, whose 
representatives are shown in Fig. [T] The mass-dependence of 



,(8) 



Af\mf_,/me,mf^/mr). 



group 


Af\mf^/me) 


A^**(mp/mT) 


AS) 
^3 


1(a) 


7.74547 (42) 


0.000032 (0) 


0.003209 (0) 


1(b) 


7.58201 (71) 


0.000252 (0) 


0.002611 (0) 


1(c) 


1.624307 (40) 


0.000737 (0) 


0.001807 (0) 


1(d) 


-0.22982 (37) 


0.000368 (0) 





11(a) 


-2.77888 (38) 


-0.007329 (1) 





11(b) 


-4.55277 (30) 


-0.002036 (0) 


-0.009008 (1) 


11(c) 


-9.34180 (83) 


-0.005246 (1) 


-0.019642 (2) 


III 


10.7934 (27) 


0.04504 (14) 





IV(a) 


123.78551 (44) 


0.038513 (11) 


0.083739 (36) 


IV(b) 


-0.4170 (37) 


0.006106 (31) 





IV(c) 


2.9072 (44) 


-0.01823 (11) 





IV(d) 


-4.43243 (58) 


-0.015868 (37) 
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FIG. 2. Self-energy-like diagrams representing 32 gauge- 
invariant subsets contributing to the lepton g — 2 at the tenth 
order. Solid lines represent lepton lines propagating in a weak 
magnetic field. 

eighth-order QED contribution a^^^: 

a|f) = -1.9106 (20) + 132.685 2 (60) 
+ 0.042 34 (12) + 0.062 72 (4) 
= 130.879 6 (63). (10) 



Over the period of more than nine years we have nu- 
merically evaluated all 32 gauge-invariant sets of dia- 

,(10) 



■22. imi 



whose rep- 
The results 



grams that contribute to 
resentative diagrams are shown in Fig. [2l 
for mass-dependent terms are summarized in Table |lll 
Some simple diagrams were evaluated analytically or in 
the asymptotic expansion in m^/me [4ll - l45l |. The results 
are consistent with our numerical ones. 

From the data listed in this Table and the value of 
A[^°^ from Ref. [13, we obtain the complete tenth-order 
result: 



a 



(,1") = 9.168 (571) + 742.18 (87) 
753.29 (1.04). 



0.068 (5) 



2.011 (10) 
(11) 



The uncertainty 1.04 is attributed entirely to the statis- 
tical fluctuation in the Monte-Carlo integration of Feyn- 
man amplitudes by VEGAS This is 20 times more 
precise than the previous estimate, 663 ( 20), obtained 
in the leading- logarithmic approximation |22| . This is 
mainly because we had underestimated the magnitude of 
the contribution of the Set Ill(a). Note also that ((TT|) is 
about 4.5 s.d. larger than the leading-log estimate. The 



numerical values of {a/ir) 



(«)„(2« 



for n = 1, 2, 



, 5 are 



summarized in Table [ 

In order to evaluate a^(QED) using ([7]), a precise value 
of a is needed. At present, the best non-QED a is the 
one obtained from the measurement of /i/mRb (47j . com- 
bined with the very precisely known Rydberg constant 
and TTiRb/we (ssj : 

a-i(Rb) = 137.035 999 049 (90) [0.66 ppb]. (12) 
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TABLE II. Tenth-order mass-dependent contribution to the 
muon g — 2 from 31 gauge-invariant subsets shown in Fig. [21 
The mass-dependence of ^3^°' is Ag^"^ (mp/me, m^/mr). 
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181.1285 (51) 


0.096 522 (93) 


0, 


,194 76 (26) 


VI(c) 


-36.58 (12) 


-0.2601 (28) 


-0, 


.5018 (89) 


VI(d) 


-7.92 (60) 


0.0818 (17) 







VI(e) 


-4.32 (14) 


-0.035 94 (32) 


-0, 


,1122 (24) 


vi(f) 


-38.16 (15) 


0.043 47 (85) 


0, 


,0659 (31) 


VI(g) 


6.96 (48) 


-0.044 51 (96) 







VI(h) 


-8.55 (23) 


0.004 85 (46) 







VI(i) 


-27.34 (12) 


-0.003 45 (33) 


-0, 


,0027 (11) 


vi(j) 


-25.505 (20) 


-0.011 49 (33) 


-0, 


,016 03 (58) 


VI(k) 


97.123 (62) 


0.002 17 (16) 








Actually, we have a more precise value of a which is de- 
rived from the measurement (isl . |49| and theory of the 
electron 5 - 2 js^]: 



\a^) 137.035 999 1736 (68)(46)(26)(331) 
[0.25 ppb] , 



(13) 



where the first three uncertainties are due to the eighth- 
order term, tenth-order term, and the hadronic and elec- 



TABLE III. Contributions to muon g — 2 from QED pertur- 
bation term a<?"^(a/7r)" x lO". They are evaluated with two 
values of the fine-structure constant determined by the Rb 
experiment and by the electron g — 2 (oe). 



order 



with Q~^(Rb) 



with OL ^(fle) 



2 116 140 973.318 (77) 116 140 973.212 (30) 
4 413 217.6291 (90) 413 217.6284 (89) 

6 30 141.902 48 (41) 30 141.902 39 (40) 

8 381.008 (19) 381.008 (19) 

10 5.0938 (70) 5.0938 (70) 

a^(QED) 116 584 718.951 (80) 116 584 718.845 (37) 



troweak terms, involved in the evaluation of a^- The 
fourth uncertainty comes from the measurement of a^- 
At present the difference between p2|) and (fT3|l is much 
smaller than the current uncertainty in the measurement 
of so that one may use either one of these two. How- 
ever, some caution must be exercised to employ Q;~^(ae) 
to calculate a^, when more accurate experiment of 
becomes available, because theoretical calculation of Og 
is strongly correlated with that of . 

Substituting ([9]), dlO]), and ((II]) in Eq. ^ and using 
P^ . we obtain 



,(QED,Rb) = 116 584 718 951 (9)(19)(7)(77) x 10" 



-14 



(14) 

where the uncertainties are from the lepton mass ra- 
tios, the eighth-order term, the tenth-order term, and 
the value of a. in (fT2l) . respectively. If we use the value 
of a in (IT3I) instead, we get 



a^(QED, fle) = 116 584 718 845 (9)(19)(7)(30) x 10"". 

(15) 

Note that the uncertainties of the lepton mass ratios, the 

eighth-order term, the tenth-order terms, and 0.(0,^) are 

improved by factors 1.7, 1.3, 20, and 1.5, respectively, 

compared with (QED, Oe) given in Eq. (99) of Ref. |50|. 

The difference between ([T4| and ([Tsl) is less than 1.2 x 

10~^^ so that we may use either one as far as comparison 

with the current experimental data is concerned. 

In view of the rather large value of fiS^^^ (TrL^J_|me) 

(12) 

one might wonder how large {m^/mf,) might be. 
As a matter of fact it is not difficult to estimate its 
size. For this purpose note that the dominant contri- 
bution to A2 {m^/me) comes from the Group IV(a) and 
the dominant contribution to A^2'^\m^/me) comes from 
the Set VI(a). Both are integrals obtained by insert- 
ing several second-order vacuum-polarization loops 112 
into the virtual photon lines of the sixth-order diagram 
A2^^(m^/me; l-l) which contains a light- by-light scatter- 
ing electron loop. Analogously the leading contribution 
to the twelfth-order term will come from insertion of 
three 112 's in ^2^' ('Ti^/mg; l-l), namely. 



^2^^'("V/'™e) 



A'i\m^/me\ l-l) 



and 



A. 



(12) 



(to^/toe) X - 0.8 X 10" 



x 10 (16) 



(17) 



noting that ^2^-' (m^/me; l-l) ~ 20 and the factor 10 ac- 
counts for the possible ways of insertion of 112. Including 
the contribution of other diagrams, the size of the 12th- 
order term might be as large as I0~^^. This is larger than 
the uncertainty of the lOth-order term in ([T4| so that it 
would be desirable to obtain at least a crude evaluation 
of this term. 
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Adding (g]), (O, and (HI]), and using a from 
(|T2|) . the theoretical value of in the standard model is 
given by 



a^(SM) = 116 591 840 (59) x 10"". 
We have therefore 



a^(exp) - a^(SM) 249 (87) x 10 



-11 



(18) 



(19) 



The size of discrepancy between theory and experiment 
has not changed much, since the tenth-order QED con- 
tribution is not a significant source of theoretical uncer- 
tainties. Let us emphasize, however, that the complete 
calculation of a|/°^ enables us to concentrate on improv- 
ing the precision of the hadronic contributions. 
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